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Abstract. We present a new algorithm to compute the integral closure of a 
reduced Noetherian ring in its total ring of fractions. A modification, appli- 
cable in positive characteristic, where actually all computations are over the 
original ring, is also described. The new algorithm of this paper has been 
implemented in Singular, for localizations of afhne rings with respect to ar- 
bitrary monomial orderings. Benchmark tests show that it is in general much 
faster than any other implementation of normalization algorithms known to 
us. 



1. Introduction 

Computing the normalization of a ring is a major tool in commutative algebra, 
with applications in algebraic geometry and singularity theory. The first general 



algorithms were proposed by Stolzenberg ( 1968 ) and Seidenberg ( 1970 1975 ). How- 



ever, the tools involved, such as extensions of the ground field and addition of new 
indeterminates, make them unsuitable for most practical applications. 

In recent years several new and more practicable algorithms using Groebner bases 
have been proposed. The basic approach, continuing the line of the works men- 
tioned before, is to compute an increasing chain of rings from the original ring to 



its normalization. This is carried out in the works of |Traverso ( 1986 ) , [Vasconcelos 



(1991 19981, Brennan and Vasconcelos (2001). To our knowledge none of these 



algorithms has been implemented and it remains unclear how practical they are. 



Also de Jong (1998), Decker et al. (1999) follow this path, applying as a new in- 



et al. 


( 


1999 


) they report an 


LAR ( 


Greuel et al. 




2009b). 



gredient a criterion for normality due to Grauert and Remmert (1971). In Decker 



It became the standard algorithm for normalization 



in computer algebra systems, being now implemented also in Macaulay2 (Grayson 
and Stillman 2009) and Magma ('Bosma et al. 1997). A good review on most of 



these algorithms can be found injSwanson and Huneke (2006 Chapter 15) 



Another approach, presented in Gianni and Trager (1997), is to use Noether nor- 
malization, reduce the problem to the one dimensional case, and apply existing 
special algorithms for that case (Ford 1987 Cohen 1993). Unfortunately, we do 



not know of any implementation of these algorithms. 



A more recent approach taken in Leonard and Pellikaan (2003) and Singh and 



Swanson (2008) is to compute a decreasing chain of finitely generated modules 
over the original ring containing the normalization. Their algorithm works only 
in the case when the base field is of positive characteristic p, where they can use 
the Frobenius map. It has been implemented in Macaulay2 and Singular, and it 
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turns out to be very fast for small p. However the computation of the Frobenius 
makes it impracticable when p is large. 

There are also very efficient methods for computing the normalization in some spe- 
cial cases. For example, for toric rings, one can apply fast combinatorial techniques, 
as explained in Bruns and Koch (2001 1. 

The algorithm we propose in this paper is a general algorithm and it is based on 
de Jong (1998) and Decker et al. (1999). In their algorithm, as we mentioned 
before, they construct an increasing chain of affine rings. They enlarge the rings by 
computing the endomorphism ring of a test ideal (see below) , adding new variables 
for each module generator of the endomorphism ring and dividing out the relations 
among them. Then the algorithm is applied recursively to this new affine ring. 
Due to the increasing number of variables and relations this can produce a big slow- 
down in the performance of the algorithm already when the number of intermediate 
rings is 2 or 3. For a larger number, it usually makes the algorithm unusable, 
as the Groebner bases of the ideals of relations grow extensively. Our approach 
avoids the increasing complexity when enlarging the rings, benefiting from the 
finitely generated A-module structure of the normalization. We are able to do most 
computations over the original ring without adding new variables or relations. 

The main new results of this paper are presented in Section [3] In Section |4] we 
describe the algorithm and show, as an application, how the i5-invariant of the ring 
can be computed. Section[5]contains several benchmark examples and a comparison 
with previously known algorithms, while Section[6]is devoted to an extension of the 
algorithm to non-global monomial orderings. 



2. Basic definitions and tools 

Let A be a reduced Noetherian ring.^ The normalization A oi A is the integral 
closure of A in the total ring of fractions Q{A), which is the localization of A with 
respect to the non-zerodivisors on A. A is called normal ii A — A. 

The conductor of A in A is C = {a G (5(A) \ aA d A} = Ann^(A/A). 

Lemma 2.1. A is a finitely generated A-module if and only if C contains a non- 
zerodivisor on A. 

Proof. If p G C is a non-zerodivisor then A ^ pA C A is module-finite over A^ 
since A is Noetherian. Conversely, if A is module-finite over A then any common 
multiple of the denominators of a finite set of generators is a non-zcrodivisor on A 
contained in C. □ 

We recall the Grauert and Remmert criterion of normality. 

Proposition 2.2. Let A he a Noetherian reduced ring and J d A an ideal satisfying 
the following conditions: 

(1) J contains a non-zerodivisor on A, 

(2) J is a radical ideal, 

(3) N{A) C V{J), where 

N{A) — {P C A, prime ideal \ Ap is not normal} 
is the non-normal locus of A. 



We assume that all rings are commutative with 1 and that morphisms map 1 to 1. 
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Then A is normal if and only if A ^ Hom/i(J, J), via the canonical map which 
maps a Cz A to the multiplication by a. 

Definition 2.3. An ideal J C A satisfying properties (1) — (3) is called a test ideal 
(for the normalization) of A. A pair {J,p) with J a test ideal and p G J a non- 
zerodivisor on A is called a test pair for A. 



By Lemma 2.1 test pairs exist if and only if A is module-finite over A. We can 
choose any radical ideal J such that p £ J \fC . 

Our algorithm computes the normalization of A when a test pair for A is known. If 
A is a reduced, finitely generated fc-algebra with fc a perfect field, then C contains 
a non-zerodivisor which can be computed by using the Jacobian ideal (cf. Lemma 
4.1 and Remark 4.6). The same holds for localizations of such fc-algebras w.r.t. any 



monomial ordering. Indeed, our algorithm is slightly more general, working when- 
ever the Jacobian ideal does not vanish. 

If A is not normal, we get a proper ring extension A C IIom^(J, J) —: A^. 

If Ax is not normal, which is checked by applying Proposition |2.2| to Ai, we obtain 
a new ring A^ by that same proposition, which then has to be tested for normality, 
and so on. That is, we get a chain of inclusions of rings 

A C ^1 C ^2 C . . . 

(with Ai — . . . , t^J/Ji, /,j ideal, and natural maps V'l : ^ ^i)- 

If at some point, we get a normal ring A^q ^ then A^q = ^ by Lemma 2.6 This guar- 



antees that if A is a finitely generated A-module, the chain will become stationary 
with Afq normal, giving an algorithm to compute the normalization. 



The fact which makes the whole algorithm practicable, is the isomorphism 

IIomA( J, J) = 1/p • {pJ -A J), 



allowing us to compute Hom^(J, J) over A. This fact, not contained in de Jong 



( 1998 ) , was found by the first named author during the implementation of the algo- 



rithm in Singular and first published in Decker et al. (1999) (see also Greuel and] 



Pfister (2008 Lemma 3.6.1) and Gianni and Trager (19971 for related statements) 



We shall prove a generalization of this isomorphism in Lemma \3.1\ which will be 
needed in the new algorithm. 

The following lemma describes the A-algebra structure of 1/p • {pJ :a J)- This will 
allow us to compute the normalization of A recursively. 

Lemma 2.4. Let A be a reduced Noetherian ring, and {J,p) a test pair for A. 
Let {u{) = P,ui, . . . ,Us} be a system of generators for the A-module pJ :a J- If 
ti, . . . ,ts denote new variables, then tj i-> Uj/p, 1 < J < s, defines an isomorphism 
of A-algebras 



Al A[ti,...,t,]//i 



1 

P 



{pJ -.A J), 



where 1% is the kernel of the map tj i— >■ Uj /p from A[ti, . . . ,ts\ — ^ -{pJ '-A J)- 



See Greuel and Pfister (2008 Lemma 3.6.7) for the computational aspects of this 
lemma. 

Example 2.5. Let / — (a;^ — y^') C k[x, y] and A = k[x, y]/I. We take the test pair 
(J,p), with ,/ :— {x,y)j^ (the radical of the singular locus of A) and p :— x (see 
Algorithm nj) . Then pJ -.^ J = (a;,?/^)^ and l/p ■ {pJ -.a J) = l/x ■ {x,y'^)^ = 



Al := A[t]/Ii where Ii 
1 1-^ y'^/x. 



— y,yt — x,y — xt)A[t]- The isomorphism is given by 



The following easy lemma gives a normalization criterion for ring extensions. It 
provides a convenient way to prove correctness of our normalization algorithm, or 
any modification, because it is independent of the intermediate steps. 

Lemma 2.6. Let tp : A ^ B be a map between reduced Noetherian rings satisfying 
the following conditions: 

(1) ij: is injective, 

(2) is finite, 

(3) B is contained in Q{ip{A)). 

Then ip induces isomorphisms Q{A) — > Q{B) and A ^ B. In particular, if B is 
integrally closed, then A is isomorphic to B. 

Proof. Since A ^ B is injective, so is Q{A) ^ Q{B) and hence A ^ B. The 
isomorphism Q{A) — ?> QiB) is clear by [sj The finitcncss of ip implies that B 
(and therefore B) is integral over A. Since ^{A) C B C _B C Q{B) = Q{tp{A)), 
we conclude that B is the normalization of ^{A), which immediately implies the 
isomorphism A ^ B. 



3. Computing over the original ring 



It has already been noticed by many authors (see for example the comments pre- 
ceding Prop. 6.65 of Vasconcelos (2005)) that the chain of rings mentioned in 
last section, or similar constructions where the number of variables and relations 



increase in each step, behaves poorly in practice. (See also Remark 5.1 ) 

There has been therefore a search for algorithms carrying out most of the compu- 
tations in the original ring. In Vasconcelos ( |2000 ), the author proposes to use 

B=\J Homs(/",/"), 

n>l 

where 5 is a hypersurface ring over which A is finite and birational and / is the 
annihilator of the S-module A/S. However, as mentioned in that same paper, 
computing B is still the hard part of the algorithm and there is no indication on 
how to do it. 



In this section we show that a chain of ring as used in Decker et al. ( 1999 ) can be 



constructed doing most computations over the original ring. In this way we obtain 
an algorithm that is usually much faster in practice. 

The purpose of this section is not only to show that the computations over the 
original ring are possible. The proofs which we provide show also how these com- 
putations can be done and thus prepare the algorithms presented in the next section. 

We start with a generalization of the isomorphism from the previous section, ex- 
pressing Homyi( J, J) as an ideal quotient, to be used later. We formulate a more 
general version than needed. For a related statement see |Swanson and Huneke| 
(120061). 



Lemma 3.1. Let A he a reduced (not necessarily Noetherian) ring, Q{A) its total 
ring of fractions, and I, J two A-submodules of Q{A). Assume that I contains a 
non-zerodivisor p on A. 



(1) The map 

$ : HomA(/, J) 



1 



(pJ -QiA) I) = J ■■Q(A) I, 



p ' ' p 

is independent of the choice of p and an isomorphism of A-modules. 



(2) IfJcA then 

PJ -QiA) I ^PJ -A I- 

Proof. (1) Let g G / be another non-zerodivisor on A. Write p — pi/po and 
q = qi/qo, with pq, go non-zerodivisors contained in A and Pi,qi £ A. 

Then c :— po^o G ^ is a non-zerodivisor and cp, cq € A with cpq G /. Since 
if G Homyi(/, J) is A-hnear, we can write 

cpifiiq) = ifiicpq) = cqip{p), 

whence ip{p)/p = ip{q)/q in Q{ A), showing that $ is independent oi p. 

Moreover, for any f (z I we have 

fjp) ^ J. ^ •fjcp) ^ ^ ^ y(cp/) ^ cpy)(/) ^ ^ J 
p cp cp cp ' 

in particular •^{p) ■ f G pJ. This shows that the image ^{ip) is in 1/p- {pJ '■q(a) I)- 
It also shows that f{p) = 'i=^ V f E I f{f) = i) <^ f = Q and hence that $ is 
injective. 

To see that $ is surjective, let q G Q{A) satisfy ql C J. Denote by ruq G Hom^(/, J) 
the multiplication by q. Then $(mg) — qp/p = q showing that <i> is surjective. 

(2) During the proof of (1) we have seen that 

pJ :q(A) I = Wip) I <y5 e HomA(/, J)}. 

Hence, the claimed equality holds if and only if f{p) G A for all G Hom^(/, J), 
which is clearly true if J C A. 

Recall the chain of extension rings from last section A G Ai G A2 G . . . We have 
seen that we can compute the normalization of A by computing the normalization 



of Ai (Lemma 2.6). The next proposition explains how to obtain a test pair in Ai 



from a given test pair in A. This is the only computation to be carried out in Ai. 

Proposition 3.2. Let A he a reduced Noetherian ring, A' = A[ti, . . . ,ts]/I' a finite 
extension ring, with natural inclusion ip : A ^ A' . If (J,p) is a test pair for A then 
setting J' = (■(/;( J))^', (J' ,ip{p)) is a test pair for A' . 

Proof. Let C be the conductor of A in Q{A) and C the conductor of A' in Q{A'). 
We know that N{A') = V{C'), N{A) ^ V{C) and ^(C) C C". Therefore V{C') C 
V(C), which proves that N{A') C N{^{A)) since iPiA) = A. We have N{A) C V{J) 
by definition of J, and hence N{A') C V{ip{J)). Now V'(p) G J' is a non-zerodivisor 
on A' and (J','0(p)) is a test pair for A'. 



Example 3.3. Recall Example 2.5 We started with A — k[x , y] / {x"^ — y^) and 
test pair {J,p) = {{x,y),x) and obtained Ai := A[t]/Ii = 1/di ■ Ui where Ii = 
{i^ -y,yt- x,y'^ - xt), di = a; and Ui = {x,y'^)A- 



We now compute Ji = ^ (V'i(J))ai = ^y{x,y)A^ = {x,y,t)Ai = (Oai (since = y 
and = a; in Ai). Therefore ((t). a;) is a test pair for Ai. 

For the remainder of this section, let i? be a Noetherian ring, / C i? a radical ideal 
and A = R/I. 

We are mainly interested in i? = k[xi, . . . , x„] with k a field (which the reader may 
assume in the following), or R = k[xi,X2, . . ■ ,Xn]> with > an arbitrary monomial 
ordering. However the proposed method works quite generally, whenever a test pair 
is known. 
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In the new algorithm, we will compute ideals Ui,U2, ■ ■ ■ , Un of A and non-zerodivi- 
sors di ^ Ui, I < i < N , on A such that 

A c ^Ui c C • • • C -I-Un = A. 

From the construction we know that 1/di ■ Ui is a finitely generated i?-algebra and 
hence there is a surjection 

1 

d,. 



R, R[ti,t2, ■ ■ . -yU„ tj Uj, 



where {di , ui , . . . , Ms . } is a set of i?-module generators of Ui . If li denotes the kernel 
of this map, we get a ring map 

LP, : A, := R,/h ^ ]rU, C Q{A). 
(Note that the definition of li is now slightly different from the one given in|2.4[) 



Example 3.4. Carrying on with Example 3.3 we compute (pi{Ji) = ipi{{t)) 



Note that (pi{t) ~ y'^/x. However the A-module (y^ Ix)a 5 'Pi{{t))- For example, 
we have seen that y e {t)Ai and clearly (pi{y) — yx/x, but yx ^ {y"^) a- 

This shows that in order to obtain A- module generators of (pi{Ji) it is not enough 
to compute the images of the generators of J^. In Algorithm [2] we will show how to 
compute the generators. In this example, it turns out that ipi{{t)) — {yx/x,y'^/x) 
as A-module. 

Once we have computed a test pair { Ji,pi) in Ai, the next step is to compute the 
quotient pJi :Ai Ji- The following theorem shows that this computation can be 
carried out in the original ring A. 

Theorem 3.5. Let A = R/I , A' = A[ti, . . . ,ts]/I' a finite ring extension and maps 
^jJ : A ^ A' , (p : A' ^ Q{A). Let {J,p) be a test pair for A and {J' ,p') a test pair 

for A' , with p' — ip{p). Let U,H be ideals of A and d ^ A such that (p(A') = -U 

and ^piJ') = —H. Then 
d 

ip'J') -.A' J' = ^{dpH -.A H). 
Proof. The proof is an easy consequence of Lemma |3.1[ Omitting ip and ijj in the 



following notations and applying Lemma 3.1 to p € J C ^ we get 

[p'J') -.A' J' = {p'J') -.QiA) J' = pH :q(a) H, 
since Q{A') = Q{A) and J' ^1/d-H. 

On the other hand, we can apply Lemma [3?T] to dp E H C A and get 



^ {dpH :a H) = -{dpH :q^a) H) ^ pH :q^a) H. 



d 



Using Theorem |3.5| together with the previous results, once we have computed an 
intermediate ring Ai, we can compute Ai^i, the next ring in the chain. If = ^i+i. 



we have finished and Ai is the normalization of the original ring A, by Lemma 2.6 
If not, we proceed by induction to compute the normalization. 

We continue with the above example. 

6 



Example 3.6. We have p ^ di = x and Hi = {xy,y^)A- We compute dipHi :j 
Hi =x2(a;y,?/2) {xy,y^) = {x^.xy'^). 

Then 



HoniAj(Ji, Ji) ^ ^{x'^,xy'^) 



1 



This is equal to Ai. Therefore, the ring Ai was already normal, and hence equal 
to the normalization of A. 

Modification 3.7. We have seen that the only computation performed in Ai is the 
radical of ipi{J). However, when the characteristic of the base field is q > it is 
possible to compute also this radical over the original ring. For this, we use the 



Frobenius map, as described in Matsumoto (2001). 



Let G ~ ipi{J) C Ai. By definition, 

J, = = {/ G A, I /" e G for some m e N}. 
Mapping to Q{A), we obtain 



Then 



d,G,-{/eC/, |/''ed,V,(G)}. 

Now d^Lpi{G) is an ideal of A and diGq is the so-called q-th root of d^(pi{G). This 
ideal can be computed over A using the Frobenius map (cf. Matsumoto (2001)). 

By iteratively computing the q-th root of the output, until no new polynomials are 
added, we obtain Lpi{Ji) as desired. 





f/d^ 




{f/d^) 


where G„j := < 






{f/d) 



Computing the radical in this way, we get another algorithm (in positive charac- 
teristic) which is similar to the one proposed in Singh and Swanson (2008). In 
their algorithm they start with the inclusion A C ^A, where c is an element of the 
conductor and compute a decreasing chain of A-modules 



-A 



\K^\u'i. 



U'^ = A. 



In our algorithm we compute an increasing chain 



Ac ^Ui c 
di 



C —Un 
dN 



A. 



The most difficult computational task for both algorithms is the Frobenius map. 
However, in our algorithm we start with a small denominator di and therefore the 
computations might be in some cases easier. This modification has not yet been 
tested. 



4. Algorithms and application 

We describe the algorithm in general terms. Since we compute an increasing se- 
quence of subrings of the integral closure the algorithm terminates, for a Noet herian 
ring A, if and only if ^ is a finitely generated ^-module. By Lemma 2.1 this is 
equivalent to the existence of a test pair. We now deal with the problem of con- 
structing an initial test pair. 
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Lemma 4.1. Let k be a perfect field, and A = k[xi,X2, ■ ■ ■ ,Xn]/I with I ~ 
if ii /2, ■ • • , /t) o, reduced equidimensional ring of dimension r. Let M be the Jaco- 
bian ideal of L, that is, the ideal in A generated by the images of the {n — r)x (n — r)- 
minors of the Jacobian matrix {dfi/dxj)i j. Then M is contained in the conductor 
of A and contains a non-zerodivisor on A. 

Proof. Let / = Pi n P2 H . . . n Ps with Pi,P2, . . . ,Ps the minimal associated primes 
of /. Since A is equidimensional, dim(A) = height(Pi) — r for 1 < i < s. Hence, the 
image of M in Ai = k[xi, X2, . . . , Xn]/Pi is contained in the Jacobian ideal Mi of Pi. 



By the Lipman-Sathaye theorem (cf. Swanson and Hunekc (2006) and Singh and 



Swanson (2008 Remark 1.5)) and hence M is contained in the conductor of Ai. 
Since A = Ai (B A2 (B ■ ■ ■ ® Ag, M is then also in the conductor of A. Moreover, the 
image of M in Ai is not zero since Ai is reduced. This follows from the Jacobian 



criterion and by Serre's condition for reducedness (cf. [Greuel and Pfister (2008 



Section 5.7)). As a consequence, AI is not contained in the union of the minimal 
associated primes of A and hence contains a non-zerodivisor on A. 

Note that both, the Lipman-Sathaye theorem and the Jacobian criterion, require k 
to be perfect. 

The ideal J := ^/M from last lemma can be used as an initial test ideal. To 
construct a test pair, we need to find in addition a non-zerodivisor of A in J. 
An element p G A is a non-zerodivisor if and only if :a (p) = 0, hence the 
non-zerodivisor test is effective. However, it is not sufficient to apply the test to 
the generators of J. (E.g., if / = (xy), the polynomials x,y generate J and are 
zerodivisors on A, but x + y is not.) Since we cannot test all elements of J there 
seems to be a problem to find a test pair if / is not prime. We address this problem 
as well as the perfectness and the equidimensionality assumptions in Remark |4.6[ 

We first describe in Algorithm [l] how to compute the initial test pair (J,p) in A, 
assuming that we are able to find a non-zerodivisor. 

Remark 4.2. Only for this step we need the assumption that R = k[xi, X2, ■ ■ ■ , Xn]> 
with k perfect and that / is equidimensional. All further steps do not require this 
assumption. 

If, by whatever means, an initial test pair {J,p) for A is known, we can start 
with the computation of Ui and then all further steps are correct, and the loop 
terminates with the computation of A. Hence, for any reduced ring A = R/I with 
given test pair {J,p), the algorithm is effective when Grobner bases, ideal quotients, 
and radicals can be computed in rings of the form R[ti, t2, . ■ . ,ts]. 



Algorithm 1 Initial test pair {J,p) 

Input: / C i?, an equidimensional radical ideal, with R = k[xi, X2, ■ ■ ■ ,a;„]> and 

k a perfect field. 
Output: {J,p) a test pair for A := R/I. 

r := dini(/) 

M' := the Jacobian ideal of /, i.e., the ideal in R generated by the 

(n — r) X (n — r)-minors of the Jacobian matrix of / 
M := the image of M' in A 
J ^/M cA 

choose p J such that p is a non-zerodivisor on A 
return {J,p) 



We now explain how to perform some auxiliary tasks, that will be needed in the 
main algorithm. 
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We have seen in the previous section that if v4 = R/I and A' ~ R[ti, . . . ,tn]/I' a 
finite extension ring with I d T , then there exist a non-zerodivisor d e yl, an ideal 
U C A and a map ip : A' ^ 1/d ■ U such that A' ^> 1/d-U. For computations, we 
need to know how to move from one representation to the other. 

Remark 4.3. If we know d and generators {d,ui, . . . ,Us} of U, we can exphcitly 
compute ip{q) for any q e A'. Let g G i?' be a representative, and substitute all 
the variables tj in q by the corresponding fraction Uj/d. This results in an element 
f/d'^G Q{A) for some f G A and e G Z>o. Now we need to find f G A such that 
//d^ ^ f/fi in Q(^)^ .^^hich is equivalent to / = fd"-'^ + g in R, with g e I. We 
can find /' by solving the (extended) ideal membership problem / € / + {d'^~^) 



in R, e.g. by using the Singular command lift, cf. Greuel and Pfister (2008 
Example 1.8.2). 

We will need also to compute ^-module generators of ideals J' C A' given by 
generators in A' . It is clear that for any such J' there exist an ideal H C A such 
that <p{J') = 1/d ■ H. So the problem is equivalent to finding elements hi, . . . ,hi 
in A that generate H as an A-ideal. In Algorithm [2] we explain how do it. 

Algorithm 2 yl-module generators 

Input: A = R/I, with R = k[xi, . . . , Xn] and I C R ideal; A' ~ R'/I' a ring 
extension of A, with R' = R[ti, . . . ,ts] and /' C R' an ideal; d G A' a non- 
zerodivisor and U' = {uq — d,ui, . . . ,Us) a such that A' = ^JJ, with map ip : 

A' A if/; J' = (A, . . . , an ideal of A'. 

Output: H = {hi, . . . , hi) A such that (p{J') = 1/d - H 
for j — 1, . . . ,m do 



compute hj such that Pi{fj) — hj/d (cf. Remark 4.3) 
end for 

set S = {hi, ... , hrn} 

for j = 1, . . . ,m; k — 1, . . . , s do 

compute hj^k G A such that hj_k/d = Uk/d ■ hj/d in Q{A) (again by Remark 



4.3) 



if hj^k ^ {S)a then 

S = SU{hj,k} 
end if 
end for 

return H := {S) 



Lemma 4.4. Let A — R/I, with R = k[xi, . . . , x.^] and I C R ideal; A' — R'/I' 
a ring extension of A, with R' = R[ti, . . . ,ts\ and I' C R' an ideal; d £ A' a non- 
zerodivisor and U' = {uq ^ d,ui, . . . , Us)a such that A' = , with map tp : A' 
^U; J' = {fi, . . . , fm) A' , an ideal of A' . The output ideal H = {hi, . . . ,hi) a of 
Algorithm^ satisfies ip{J') = l/d ■ H. 

Proof. This follows since the ^-module (1 = uo/d,ui/di, . . . ,Us/d)j^ = piA') and 
the A'-module {hi/d, h^/d, . . . , hm/d)^, = (p{J') {hi, . . . ,hm as in the algorithm). 
Therefore the products Uk/d ■ hj/d, < k < s,l < j < m, generate ip{J') as 
A-module. Hence {hj \ 1 < j < 1} generates H as ^-module, or equivalently as 
A-ideal. 



Example 4.5. We apply the algorithm to compute the A-module generators of 
Lpi{Ji) from Example 3.4 Recall that Ji = {t)Ai, Ui = {x,y^)A and d — x. 
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We start with hi — filt) — jx and S = {hi\. In the first step, we compute 
x/x ■ jx — /x, therefore — ■ Since G (2/^)1 we do not do anything. 
In the second step we compute jx ■ jx = y"^ jx^ — x^yjx^ = xy/x, therefore 
^1,1 = xy- Since xy ^ (y^), we add it to S. We finish with H = {xy.y"^), as 
mentioned in Example |3.4[ 



We are now ready to present in Algorithm [3] the main algorithm to compute the 
normalization. 



Termination follows from Lemma 2.1 and the discussion after Definition 2.3 
rectness follows from Lemma [SHI 



cor- 



Algorithm 3 Normalization of R/ 1 



Input: / C i?, an equidimensional radical ideal 

Output: generators of an ideal U C R, and d E R such that A = -U C Q{A), 

with A R/I. 
compute { J,p), an initial test ideal 
C/i {pJ ■.aJ)cA 
di p 

if {di) = Ui then 

return 
end if 
i := 1 
loop 

write Ui = (^di,Ui \u2 \ ■ ■ ■ '"s*"*^^ 

set Ri R[ti, . . . ,1^], -ki : Ri ~^ j^Ui (1 ^A the map tj i-> u'^^ jd^ 



4.41 



li :— ker(7ri) (cf. Lemma 2.4 1 
set Ai = Ri/Ii 

A ■■= \fh{J) C A,, with tPi-. A^ Ai 
compute {/i, . . . , /fc} C A such that (/i, /2, ■ ■ ■ , fk) a = di(pi{Ji), 

with ipi : Ai 'jrUi (cf. Lemma 
compute generators of Ui^i := {pdiHi 
if diUi C Ui+i then 

return {Ui,di) 
end if 

i := i + 1 
end loop 



Remark 4.6. Let us comment on some variations and generalizations of Algorithm 
|3] For this let k be any field, R = X2, . . . , x„]>, and I d R a radical ideal. 
(1) If / is not (or not known to be) equidimensional we can start with an algorithm 
to compute the minimal associated primes (cf . [Greuel and Pfister ( 2008 Algorithm 



4.3.4, Algorithm 4.4.3)) or the equidimensional parts (cf. Greuel and Pfister (2008 



Algorithm 4.4.9)) of /, where the latter is often faster. The corresponding ideals 

Ji, J2, . . . , /r are equidimensional and we have R/I = R/Ii (B R/h ® ■ ■ ■ ® R/Ir- 
Hence the problem is reduced to the case of / being prime or equidimensional. 
(2) Now let / be equidimensional and M the Jacobian ideal. Since regular rings 
are normal, it follows from the Jacobian criterion that N{R/I) C V{M). Let us 
assume that M 7^ and choose p £ M \ {0}. 

a) If Ii := I :r (p) C / then p is a non-zerodivisor on A and J = ^/M is a test 
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ideal. This is always the case if / is prime. 

b) li Ii <^ I we co mpute I2 := I :r Ii and obtain a splitting I = Ii O I2 (cf. Greuel 
and Pfister (2008 Lemma 1 



.14(3))) and R/I = R/h ® R/h- Hence we can con- 
tinue with the ideals Ii and I2 separately which have both fewer minimal associated 
primes than /. Consequently, after finitely many splittings, the corresponding ideal 
is prime or we have found a non-zerodivisor. This provides us with test ideals as 
in case a). 

(3) The above arguments show that (even if k is not perfect) Algorithm [3] works 
for prime ideals if and only if the Jacobian ideal M is not zero. This is always the 
case for k perfect. However, if k is not perfect, M = may occur. For example, 
consider k = {Z/q){t) with q a prime number, and / — {x^ + y'^ + 1) C k[x,y]. For 
a method to compute a non-zero element in the conductor oi R/I ii I is prime and 
if Q{R/I) is separable over k, see Swanson and Huneke (2006 Exercise 12.12). 



4.1. The (5-invariant. As an application of the normalization algorithm we show 
how to compute the (5-invariant oi A = k[xi^X2, ■ ■ ■ ,Xn\>/I, a reduced Noetherian 
fc-algebra, 

S{A) :=dimfe(A/v4). 

S{A) may be infinite but it is finite if the algebraic variety V{I) defined by / 
has isolated non- normal points, e.g. for reduced curves, i.e. dim(yl) = 1. In this 
case, 6 is important as it is the difference between the arithmetic and the geo- 
metric genus of a curve. Moreover, the J-invariant is one of the most important 



numerical invariants for curve singularities (cf. Campillo et al. (20071), that is, for 
1-dimensional complete local rings A. The extension of our algorithm to non- 
global orderings in Section [6] has the immediate consequence that it allows to 
compute S for affine rings as well as for local rings of singularities, noting that 
S{k[xi, X2, . . . , Xn]>/I) — 5{k[[xi, X2-, ■ ■ ■ , x„]]//) if > is a local ordering. 

Lemma 4.7. Let R he a reduced Noetherian ring, I <Z R he a radical ideal, and 
I = Pi n ■ ■ ■ n Pr its prime decomposition. Write I = Ii H ■ ■ ■ H Is, where Ii = 
CljeN "^^"^ {^ii • • ■ I -^s} is a partition of {1, . . . , r}. Let Ui, di he the output of 
the normalization algorithm for Ai = R/Ii- Then 

(1) 5{A{) = dimk{U,/d,U,), l<i<s, 

s s—1 

(2) (5(i?//) = ^5(A,)+ ^dimfe(i?/(/ + /«)), «;/iere/W =/,+in 



In particular 5 (R/I) < 00 iff every summand on the right hand side of is finite. 



Proof. This follows by induction on s, and by repeatedly applying the following 
sequence of inclusions for 5 = 2, i.e. I = Ii D I2, 



R/I R/h ® R/I2 ^ R/h ® R/h = R/h 
and the exact sequence 







R/I ^ R/h © R/h R/ih + h) 



0. 



Note that dim^ i?/(/, -f /(*)) can be computed from a standard basis of h + J^'^ 
and (lim.k{Ui/ diUi) from a standard basis of a presentation matrix of U/diUi via 
modulo (cf. Greuel and Pfister (2008 Singular Example 2.1.26)). An algorithm 
to compute 6 is also implemented in Singular (Greuel et al.l 2009b). 
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Table 1. Timings 



iNO. 


char 


normal data 


seconds 


non-zerodivisor 


steps 


normal 


normalP 


normalC 


1 





y 


7 





_ 


72 


1 


2 


y 


7 











1 


5 


y 


7 


1 


73 





1 


11 


X - 2y 


7 


1 


12 


* 


1 


32003 


y 


7 







1 


2 





y 


7 


1 


- 




2 


3 


y 


8 








3 


2 


13 


y 


7 







10 


2 


32003 


y 


7 







10 


3 





y 


6 


2 


_ 




3 


2 


y 


13 


1 







3 


5 


y 


6 


1 


7 




3 


11 


X + 4y 


6 


1 






3 


32003 


y 


6 


1 






4 





2x'^y -y^ +y 


1 





- 





4 


5 


x^y + 2y^ - 2y 


1 





3 





4 


11 


x'^y + 5?/^ - by 


1 










4 


32003 


.T^y+16001y^-16001j/ 


1 





* 





5 





y 


1 










5 


5 


x'^y + xy 


3 


1 




* 


5 


11 


y 


1 











5 


32003 


y 


1 


1 


* 





6 


2 


V 


2 


6 


24 


172 


7 





y 


6 


12 




582 


7 


2 


y 


6 


11 





35 


7 


5 


y 


6 


12 


3 


358 


7 


11 


y 


6 


11 


43 


503 


7 


32003 


y 


6 


11 


* 


617 



5. Examples and comparisons 



In Table [T] we see a comparison of the implementations in Singular of the new al- 
gorithm normal and other existing algorithms. normalC is an implementation based 



on the algorithm Decker et al. (1999) (see also Greuel and Pfister (2008 Section 



3.6)) and normalP is an implementation of the algorithm of Leonard and Pellikaan 
(2003), Singh and Swanson (2008) for positive characteristic. All these implementa- 



tions are now available in the Singular library normal . lib ( Greuel et al. 2009a ). 



Computations were performed on a compute server running a 1.60GHz Dual AMD 
Opteron 242 with 8GB ram. 

* indicates that the algorithm had not finished after 20 minutes, 

- indicates that the algorithm is not applicable (i.e., using normalP in characteristic 

0). 

We try several examples over the fields A; = Q and k = S {2, 3, 5, 11, 13, 32003}, 
when the ideal is prime in the corresponding ring. We see that the new algorithm 
is extremely fast compared to the other algorithms. Only the algorithm normalP 
is sometimes faster for very small characteristic. 
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In columns 3 and 4 we give additional information on how the new algorithm 
works. The column "non-zerodivisor" indicates which non-zerodivisor is chosen. 
The column "steps" indicates how many loop steps are needed to compute the 
normalization. We see that our new algorithm performs well compared to the 
classic algorithm especially when the number of steps needed is large. 

We use the following examples: 

• /i {{x - y)x{y + x^f - y^{x^ + xy - y^)) C k[x, y], 

• l2 = (55x« + 662/2 a;9 + 837a;2/ _ 75y^x^ - 70y^ - 97y^x'^) C k[x, y], 

• I3 = {y^ + y^x + y^ +y^ + y^x + y'^x'^ + y^x^ + yx^ + x'^) C k[x, y], 

• /4 = +y2 _ 1)3 + 27x2?;2) c k[x,y\, 

• h = + x^y^ - x^y'^ - x'^y^ + 2y^'^ ~ x^ + 2x*^y'^ + x'^y'^ ~ x^y^ - ?/ + 
2x^ - x^y"^ + x'^y^ + 2.x'' + 2x2y'^ - y'^ - x^ + y'^ - I) C k[x,y], 

• Iq = {z^ + zyx + y'^x^ + y^x^ , uyx + z'^ , uz + z + y^x + yx^, + u + zy + 
zx, v'^ +VUX + uz^ +vzyx + vzx + uz^ + uz'^y + z"^ + z'^yx'^) C k[x,y,z,u,v\. 

• Ij — {x"^ + zw,y^ + xwt,xw^ + z^t + ywt^,y^w'^ — xy^z^t — w^t^) C 
k[x,y,z,w,t]. 

Remark 5.1. As mentioned before, the main drawback of the algorithm [Decker] 



et al. ( 1999 ) is the increasing complexity of the new rings that are constructed. A 
direct implementation of the algorithm turns out to be so slow that it does not even 
finish for most of the examples analyzed in this paper (after 1 hour) . For example, 
in the second example (12) over Z3, the fifth ring constructed in the chain has 12 
variables and 76 generators for the ideal of relations. The sixth ring could not be 
computed using this direct approach. 

A partial solution to this problem, used in implementations, is to eliminate as far 
as possible redundant variables, that is, variables than can be expressed in term of 
the others through the relations in the ring. This is what is done in normalC, and 
it is sometimes a good improvement. However detecting the redundant variables 
becomes more and more difficult as the relations get more and more complex, 
adding a new expensive task to the computation, that does not always succeeds in 
detecting all the relations. 

The algorithm proposed in this paper avoids this problem in a natural way. 

We have also compared our implementation with the normalization procedures in 



Macaulay2 (they use the algorithms Decker et al. ( 1999 ) and [Singh and Swanson 



(2008)) and in Magma (they say that they use Decker et al.| (1999) for the gen- 
eral case; however it seems to work only in characteristic and the code is not 
accessible). Our new algorithm is always faster and succeeds where the other im- 
plementations do not finish. Wc do not know implementations in other computer 
algebra systems. 

6. Extension to non-global orderings 

In this section, let > be any monomial ordering on the set Mon(xi, . . . , a;„) of 
monomials in x = (xi, . . . , a;„). That is, > is a total ordering which satisfies 

Va, P,-feZl > ^ x°'+'' > x^+'-', 

but we do not require that > is a well ordering. The main reference for this section 



is Greuel and Pfister (2008) where the theory of standard basis for such monomial 
orderings was developed. 

We consider the multiplicatively closed set 

S> := {u e k[x] \ {0} I LM(u) = 1}, 
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where LM denotes the leading monomial. The localization of k[x] w.r.t. S'> is 
denoted as 
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k[x]> Sy^k[x] = |; 



f,ue k[x],LM{u) = 1 



It is shown in Greuel and Pfister (2008 Section 1.5) that fc[a;]> is a regular Noe- 
therian ring satisfying 

k[x] C k[x]> C k[x]^,,), 

where k[x] ^^.^ denotes the localization of k[x] w.r.t. the maximal ideal {xi , . . . , Xn) = 
(x) . Note that 

• = fc[x] 4^ > is global (i.e. Xi > 1, i = 1, . . . , n), and 

• fc[a;]> = fc[x]^j.^ > is local (i.e. Xi < I, i — 1, . . . ,n). 

In applications, in particular in connection with elimination in local rings, we need 
also mixed orderings, where some of the variables are greater than and others 
smaller than 1. An important case is the product ordering > — (>i,>2) on 
Mon(a;i, . . . , a;„, t/i, . . . , j/m) where >i is global on Mon(2;i, X2, ■ ■ ■ , Xn) and >2 is 
arbitrary on Mon(yi, ?/2, ■ • • , Vm)- Then 

k[x,y]> = {k[y]>^)[x] = fc[y]>., k[x], 

(cf. Greuel and Pfister ( |2008[ Examples 1.5.3)), which will be used in the extension 
of our algorithm to non-global orderings. 

We now show that for any monomial ordering > and any radical ideal / C A;[a::]>, 
the normalization of the ring k[x]y/I is a finitely generated fc[a;]> //-module and 
how to extend Algorithm |3] from Section |4] to this general situation. 

For any ideal / C k[x]y we have / = I'k[x]y, with I' = I D k[x]. Let {k[x]/I')-^ 
(resp. {k[x]/I')-^) denote the localization w.r.t. the image of S'> in k[x]/I' (resp. in 
k[x]/r). We have fc[a;]>// ^ (fc [a;] //')>. 

Lemma 6.1. With the above notations, we have an isomorphism 



k[x]>/I^{k[x]/I') 



of k[x]^- algebras. In particular, k[x]^/I is a finitely generated k[x]y/ 1 -module. 



Moreover, let k[x]/I' = k[x, t\/H as k[x]-algebras with new variables t — (ti, . . . ,ts) 
and H an ideal in k[x,t]. Then 



k[x]^/I = {k[x]y)[t]/ H{k[x] 



Proof. The first statement follows immediately from the well-known fact that lo- 
calization commutes with normalization. Since k[x]/I' is module-finite over k[x]/I' 
the same holds for the localization {k[x]/I')^ over {k[x]/I')y. The last statement 
follows since the image of S'> in fc[x,i] localizes k[x,t] only w.r.t. the x variables. 

Remark 6.2. Let /i, /2, ■ ■ ■ , fs G k[x] generate / — (/i, /2, ■ ■ • , fs}k[x]y and let /' 
denote the ideal generated by /i, /2, . . . , /s in A: [a;]. We can compute /c[a;]>// in two 
different ways. 

The first method is to compute a test ideal J and Ilomk[x]^ /i{J, J) in the same 
manner as described in the previous sections, just w.r.t. the ordering >, i.e. in 
fc[x]>. When adding new variables ti (corresponding to fc[a;]>-module generators 
of Homfe[j.]^//(J, J)) we define on k[t,x] a block ordering (>i, >) with >i a global 
ordering on the (first) t-block (i.e. ti > I for all i and ti > Xj for all and > the 
given ordering on the (second) a:-block. Then we continue with this new ring and 
monomial ordering. 
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This algorithm is correct (by applying Lemma |2.6| to A — fc[x] > // ) and terminates 
because /c[x]>// is finitely generated over fc[a:]>// by Lemma 



6.1 



The second method is to compute the normalization of k[x]/r as in the previous 
section, with all variables greater than I. Then we map the result to k[t, 



\>i,> 

with block ordering (>i, >) as for the first method. By Lemma 6.1 both methods 
give the same result, hence the second algorithm is also correct. 

If we start with an equidimensional decomposition /' = Hi^i then of course we 
only need to compute the normalization for those ideals li for which a standard 
basis of li w.r.t. the ordering > does not contain 1. 

Example 6.3. To see the difference between both methods, let 

(2/2_a;2(x + l)2(a; + 2)) Ci?:-fc[x,2/]>, 

with > a local ordering (i.e. fc[x,y]> = k[x,y](^x,y))- Let I' — I C\ k[x,y\. In Figure 
[T]we can see the real part of the curve V(J'). This curve has two singularities, at 
the points Pi = (0,0) and P2 = (-1,0). 




Figure I. y"^ ~ x'^{x + lf{x + 2) 



We carry out the first method, setting A ~ R/I. The singular locus of / is J = 
{x,y), which is radical. This is the first test ideal. We take as non-zerodivisor 
p '■— y and compute the quotient 

Ui :== yJ -.A J = {x,y). 

Since Ui ^ (y) we go on. The ring structure of 1/y ■ Ui is Ax = a;,y]>i^>//i, 
with block ordering (>i, >) (>i any ordering) and /i = (tx* + Atx^ + btx^ + 2tx — 
y, ~ty + X, t^{x + lf{x + 2) - 1, a;'^ + 4a;'' + Sx^ + 2x'^ - y^). 

We compute Ji := ^yipi{{x,y)) = {x, y, 2t^ - 1) Ai ■ 

Mapping Ji to Q{A) using di — y as denominator, we get Ji = 1/y ■ Hi, with 
Hi := {yx, y^). (The image of 2i^ — 1 in Q{A) is (— lOa^y — 8x^y — 2x'^y)/y, which 
is already in 1/y ■ {yx,y'^).) We compute the quotient 

U2 := y'^{yx,y'^) -.a {yx,y'^) {xy,y'^). 

We see that yUi = 1/2- This means that Ai was already normal and isomorphic to 
the normalization of A, which is therefore 1/y ■ {x,y)A- 

Let us now apply the second method. We set R' :— k[x,y] and A' = R'/I'. The 
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singular locus of /' is J = {x^ + x,y), which is radical. J serves as first test ideal. 
As non-zerodivisor we choose p := y and compute the quotient 

Ui := yJ :a' J = {y, + 3a;^ + 2x). 

As Ui 7^ (y), we continue. We compute A'l, the ring structure of 1/y ■ Ui, A'^ = 
k[t, X, y]/{tx^ + tx-y, -ty + x^ + Sx"^ + 2x,t^ -x- 2, x^ + Ax^ + 5x^ + 2x^ - y^), 
and Ji = y^v?i((a;2 +a;,2/)) = ^3,2 _^ ^^^y 

Mapping Ji to Q{^A!^ using rfi = y as denominator, we obtain Ji = \/y ■ Hi, with 

Hi := {y{x'^ + x),y'^). We compute the quotient 

U2 -y^yix^ + x),y^) -.a' {y{x^+x),y^) = {y^ ,y{x^ + 3x^ + 2x)). 

Now we have yUi = U2, and thus A'l was already normal and isomorphic to the 
normalization of A'. Therefore, the normalization A equals l/y-{y, x^+3x^+2x)a = 
^/y- {y,x)A, as before. 

Remark 6.4. In the previous example, using the first method yields simpler test 
ideals and quotients. However, our experience is that in general, computations 
with non-global orderings are often slower than computations with global orderings, 
and therefore the second method should be preferred at least if the input ideal is 
prime. On the other hand the computation should be faster with the first method 
if the ideal, or its jacobian ideal, has complicated components which vanish in the 
localization. 
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